Interaction-free measurements with superconducting qubits 
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An interaction-free measurement protocol is described for a quantum circuit consisting of a su- 
perconducting qubit and a read-out Josephson junction. By measuring the state of the qubit one 
can ascertain the presence of a current pulse through the circuit at a previous time without any 
energy exchange between the qubit and the pulse. 
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We know from classical optics that when an object 
blocks one arm of an interferometer the fringes disap- 
pear. In quantum physics, this effect acquires far more 
subtle features, and, counter-intuitively, it does not orig- 
inate from the unavoidable disturbance on the photon 
trajectory or random phases due to scattering, but from 
the possibility of obtaining which-path information. In- 
deed, in the quantum treatment the interference would 
disappear if one can in principle obtain information about 
which path the light went through, even if this informa- 
tion is not actually extracted Quantum processes 
therefore result not only in changes in observables such 
as position and momentum (due to energy exchange), 
but, more importantly, in establishing new correlations 
between parts of the system. Recently, a lot of work 
has been put into harnessing the power of these correla- 
tions for performing computational tasks which are very 
difficult to implement on classical computers. Supercon- 
ducting qubits based on the Josephson effect have been 
proposed [1] as the elements of future quantum comput- 
ers, based on macroscopic quantum coherence effects in 
charge and flux devices Several species of supercon- 
ducting qubits are currently under study, for example 
charge qubits [1], phase qubits [H, flux qubits 0, and a 
mixed charge-flux version called Quantronium [7|. This 
last type has a very large decoherence time (more than 
500ns), and it will be the main focus of this paper. Be- 
sides quantum computing, fundamental research such as 
testing quantum mechanics at the macroscopic level is 
an important direction envisioned decades ago with 
progress in this direction now enjoying a firm experimen- 
tal basis. 

In this paper we propose an experiment in which a 
Quantronium device (Fig. [T]) could be used to ascertain 
the presence of a small pulse of electric current without 
any disturbance due to energy exchange with the con- 
tinuum of states outside the washboard potential well 
in which the qubit is localized. An experiment of this 
type is feasible with the current Quantronium setup, and 
it would constitute a test, at the macroscopic level, of 
a strongly nonclassical prediction of quantum mechan- 
ics. The proposal is based on the interaction-free mea- 



surement scheme proposed by Elitzur and Vaidman for 
optical Mach-Zehnder interferometers Q. This realiza- 
tion has found interesting applications in interaction-free 
imaging of objects {loi |. and more recently in quantum 
computing 

In the optical setup (upper schematic in Fig. [5]), a bal- 
anced Mach-Zehnder interferometer is constructed from 
two 50% beam splitters, two perfectly reflecting mirrors, 
and two detectors D+ and D— {+ and — are the di- 
rections corresponding respectively to paths along the 
upper and lower arms of the interferometer). In the ab- 
sence of an object in the lower arm this arrangement 
produces a destructive interference at the detector D+ 
(the state of the photon, which is initially \+ >, becomes 
(I -I- > >)/y/2 inside the interferometer and i|— > 

after the second beam splitter). If a quantum ultrasen- 
sitive object {i.e. triggered by the absorption of a sin- 
gle photon) is present, in 50% of the cases there will be 
no absorption event: the photon which has traveled in 
the upper arm of the interferometer must be then in the 
state |-f > (it did not "interact" with the object), and 
it will emerge from the second beam splitter in the state 
(1+ > >)/V2, having now a 50% chance of being 

detected by D+ |l9|]. As a result, when D+ is triggered 
we can be certain about the presence of an object in a 
region in space without having to exchange energy {e.g. 
by photon absorption) with that object. The success rate 
(the fraction of photons detected by D+) is 25%. 




FIG. 1: Quantronium circuit with flux-current compensation. 

A Quantronium device (Fig. [T|) consists of a split 
Cooper pair box of total capacitance Cs and Josephson 
energy Ej/2 per junction, operated at a gate voltage Vg 
and excited by microwave radiation coupled through the 
gate capacitance Cg. The box is connected in parallel 
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with a read-out Josephson junction Ejp of capacitance 
Co- The Hamiltonian of the circuit is 7] 



H = 



1 



{q-CgVgf-Ej COS 
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Ejo cos 7 - I(l)o'j, 
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2Cs 

where Q and 7 are the charge and the phase across the 
large junction, / is the bias current, (j}4>o is the externaUy 
apphed magnetic flux through the loop, ^0 — 7i/2e is 
the reduced flux quanta, and {q,0} is the pair of conju- 
gate variables (commutation [(j)od,q] = ih) correspond- 
ing to the split Cooper pair box. The read-out junction 
is in the large-capacitance regime, Cq » e^Ej^, where 
Q becomes a continuous operator Q = —2eid/d^. A 
two-level system (the qubit) is realized 0, 3] at the gate 
voltage Vg = e/Cg, where the charging-energy degener- 
acy of the eigenvectors |0 > and |1 > of the operator q 
[q\Q >= 0, g|l >= 2e|l >) is lifted by the tunneling term 
— Ej cos {j /2+(j)/2) cosO. The eigenvectors of this tunnel- 
ing term (the qubit levels) are denoted by |+ > (ground 
state) and |- > (excited state), |± >= (|0 > ±|1 >)/a/2, 
and in this basis the Hamiltonian Eq. ^ is of Stern- 
Gerlach type 

2e2 92 - -Ej 7 + , , 

H ^ - — -^-EjoCOS-f-I(f>Qj-—cos^—cT^. (2) 

The bias current and the flux are externally controlled 
parameters {4>{t), I{t)) which are manipulated adiabati- 
cally compared to the time-scale of the qubit frequency. 
The eigenvalue-eigenfunction problem Hil}±{'-^)\± >= 
£±'0±(7)|i > can be thus solved at every moment t, 
where, to keep the notation simple, we will specify the 
corresponding values of (0, /) each time we refer to e± or 
to the instantaneous eigenvector ip±{'j). For example, at 
I = (j) — the states ^p± (7) do not differ much from one 
another, the phase 7 is almost classical- like (7 w 0), and 
the qubit energy is e = e_ — e+ w Ej Adiabatic ex- 
cursions at nonzero values of / help differentiate between 
the states ±: a typical measurement protocol for this cir- 
cuit proceeds by raising the bias current to a value close 
to the critical current EjQcf)^^ of the large read-out junc- 
tion for a certain read-out time r^. The junction then 
can tunnel in the running-wave state [l^ with a switch- 
ing rate 0] which depends on the state |± > of the 
qubit . 

The idea of the proposed experiment is to create the 
time-equivalent of a Mach-Zehnder interaction-free ex- 
periment by inserting a bias pulse (referred to in the fol- 
lowing as " interaction pulse" ) inside a Ramsey sequence 
of n/2 microwave gate pulses. Ramsey techniques are 
common to many fields of physics; they can be inter- 
preted as two-paths interferometry, the first tt/2 pulse 
separating the paths and the second rejoining them to 
test what has happened in between. The experiment pro- 
ceeds as follows (Fig. [5]): at all times, the system is kept 



at the charge degeneracy point CgVg = e where the de- 
coherence due to gate voltage fluctuations is zero in the 
first order. A first 7r/2 pulse initializes the qubit in an 
equal- weight superposition of the states |-|- > and |— >. 
After that the bias current is increased adiabatically to 
a value I*-^-' that allows tunneling during a time interval 
Tp, then turned off back to zero. Finally, another tt/2 
pulse is applied and after that a standard switching cur- 
rent measurement sequence 0, Il3| follows. Therefore, in 
this proposal the role of the ultrasensitive quantum ob- 
ject is played by the quasi-continuum of modes outside 
the well of the washboard potential of the large junction, 
into which the phase 7 is allowed to tunnel during the 
interaction pulse. 
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FIG. 2: Mach-Zehnder interferometric setup (upper 
schematic) and the Quantronium pulse sequence (graphs be- 
low) with the axes of the flux graph slightly displaced for 
clarity. 

The duration of the whole procedure must be less than 
the decoherence time T2 ~ 500ns, and the value of i^^-* 
has to be such that the dominant switching process is 
macroscopic quantum tunneling (and not thermal acti- 
vation). It is simple to maintain the qubit at the op- 
timal charge parameter point CgVg = e, but changes 
in the bias current will displace the system along the 
phase direction We propose to compensate this dis- 
placement by adding a magnetic field pulse such that 
(f){t) = — 7m(i), where jm(t) is the phase corresponding 
to the minimum of the washboard well potential of the 
large junction at a time t during the current pulse, de- 
fined by sin7„(i) = I{t)4>o/Ejo, < 7^(t) < 7r/2. As 
a result, there will be no linear longitudinal noise com- 
ponent due to fiuctuations of the fiux in the loop {e.g. 
vortices moving in and out, fluctuations of the externally 
applied magnetic field) and the decoherence time T2 is 
expected to stay of the order of 0.5/^s. Technically, it 
is easier to apply a magnetic field pulse and then model 
the current bias pulse according to the relation above 
{I(t) = — i?jo0o ^ sin0(i)), as shown in Fig. [H This si- 
multaneous manipulation of the bias current and fiux re- 
sults in a different physics than that of the usual read-out 
pulse. In the last case, the two states of the qubit are dis- 
tinguished by the appearance of small electrical currents 
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that substract from or add to the externally imposed bias 
current. In our procedure, during the excursion toward 
the switching point, these currents are maintained to zero 
to reduce decoherence; however, the two states are still 
distinguishable by different values of the plasma oscilla- 
tion frequency and of the tunneling barrier height, which 
lead to different macroscopic quantum tunneling proba- 
bilities. Finally, one can notice that the qubit Larmor 
frequency e remains unchanged during the compensated 
interaction pulses, therefore no extra phase differences 
are introduced in the Ramsey interference pattern. 

At the plateau 7*^^^ of the interaction pulse, 
by expanding the potential energy around 7m'' — 
arcsin (^oZ-E-Jo)! we find the Josephson plasma fre- 
quency 

and the corresponding washboard potential height 
n 3 



read-out pulse for the state Eq. ([31), with the result 



(2/3)^70 
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WKB approximation 
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We start with the qubit in the ground state >; 
after the first tt/2 Ramsey pulse this state (all the 
states from now on will be written in the interaction 
picture) changes into the superposition {ip+{'j)\+ > 
+iip-{'-/))\— > /V2 [H, 0]- In the absence of an inter- 
action pulse, one can check that after the second Ram- 
sey 7r/2 pulse the state of the system is itp-{'y)\— > 
[l| (destructive interference toward the D+ detector, in 
the Mach-Zehnder setup). Therefore, during the read- 
out pulse of length Tr, the non-switch probability is 
exp(-FL''V^). 

In the case in which we have an interaction pulse, 
the status of the quantum circuit immediately at 
the end of the interaction pulse depends on whether 
the large junction has switched or not. Immedi- 
ately after the interaction pulse, if the large junc- 
tion did not switch, which happens with probabil- 
ity TV = [exp(-F^Vp) -f exp(-F^Vp)j /2, the state 



r(p) /%, ,(P) 



of the circuit is (27V)~i/2 exp(-F^Vp/2)^+(7)| + 

-fz(2AA)-i/2exp(-FL^)Tp/2)V'-(7)|- >■ 
Ramsey pulse transforms this state into 



> 



The second 
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e + - e 
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r<i''Tp/2 



^-+(7)1+ > 
)V-(7)|->- (3) 



Finally, after the read-out pulse, the probability that the 
system did not switch during the entire sequence of pulses 
is a conditional probability obtained by multiplying the 
probability J\f that the system did not switch during the 
first pulse with the nonswitch probabilities during the 



/2_e-rL^'T^/2\-g-r-'r, 



/4 



/2 ,-ri->rj2Y -r^"^ 
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We analyze now the read-out pulse sequence. To avoid 
the case in which non-switching events occur in the ab- 
sence of the interaction bias pulse, we have to set the 
corresponding probability exp(— f'^V^''^) as small as pos- 
sible: suppose we consider this spurious non-switching 
rate as satisfactory when it is of the order of 10~^. 
To determine the read-out bias pulse parameters, we 
will use directly the experimental data reported so far 
Suppose we use a bias read-out current of 1.11/iA 
(94.9% of the critical current EjQ(j)Q^); then at 14 mK 

(r) 

the switching can be estimated as F_' — 15MHz, and 
F^''^ = 2MHz. For a read-out time t^'') = 460ns we ob- 
tain exp(-FL''V('')) = 10-3 and exp(-F^''V('')) = 0.4. 

Let us now focus on the interaction pulses. For an in- 
teraction bias pulse of height I^p^cJiq/Ejo = 95.2%, we 



obtain F 



(p) 
+ 



2.22MHz and F 



(p) 



6.58MHz. Since 



exp(— F^lV^''^) is negligible, we find from Eq. ([4]) that 
the maximum success rate is obtained when the differ- 
ence between the nonswitching probabilities of the two 
states, exp(— F^Vp/2) — exp(— FffVp/2) is maximal. For 
Tp « 250ns - a value of the same order as that used 
routinely in switching current experiments - we reach 
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exp(-F^Vp/2) - exp(-F^Vp/2) /4 = 2.5%. Note 

that the pulse duration is below 500 ns; we have checked 
also that the thermal activation rate is negligible. Now, 
taking into account the exp(— F^'r''''') = 0.4 read-out 
discrimination factor estimated above, we obtain as the 
final result a success rate of 1%, much larger than the 
error due to spurious nonswitching events and easily de- 
tectable experimentally. In a typical experiment record- 
ing about 10* events per second, this corresponds to 
100 successful events every second. This success rate is 
smaller than the theoretical maximum of 25%, result- 
ing from the imperfect discrimination between the states 
1+ > and |— >. Due to this, our estimates show that the 
success rate cannot be improved by using the quantum 
Zeno effect In a Mach-Zehnder setup, this would cor- 
respond to the quantum ultrasensitive object interacting 
with different probabilities with photons propagating in 
both arms of the interferometer. 

It is instructive to consider the case in which a phase 
error A = J 6e(t)dt {6e(t) is the change of the qubit 
energy at each instant t) is introduced due to an im- 
perfect compensation (f){t) ^ —"fm(t). This will pro- 
duce a dephasing of the qubit with respect to the mi- 
crowave radiation at the second tt/2 Ramsey pulse. Im- 
mediately after the interaction pulse, the circuit is in 
the state (27V)-i/2 exp(iA/2) exp(-F^Vp/2)i/)+(7)|+ > 
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+i(2Ar)-i/2exp(-iA/2)exp(-rL^Vp/2)V'-(7)|- >, and 
after the second Ramsey pulse the state assumes the form 

The probability of not observing a switch dur- 
ing the entire sequence (with the read-out strategy 
exp(— r'^V'^'"^) « adopted before), is obtained as 

(e-i^<''W2_,^r<^V,/2^)%-r<:'.(^./4 

+e-(^^'+^'-')^^/^sin^(A/2)e-^t'^'r (5) 

One can now clearly distinguish between classical inter- 
ference (oscillations of probability as a sine function, the 
second term of the expression in square brackets above) 
and the pure interaction-free effect (the first term in 
the sum above). Classical optical interferometry allows 
the detection of a transparent object due to its index 
of refraction (resulting in different optical path lengths). 
Interaction-free detection allow us to assert the presence 
of perfectly opaque objects without any photon being 
absorbed. In general, if A is nonzero, both effects are 
present, as shown in Eq. ([5|). In the case of the parame- 
ters suggested above for the interaction pulse, a dominant 
interaction-free effect is obtained if sin^(A/2) <c; 7.5% 
(within the reach of present-day experiments 0). The 
experimental setup can be also calibrated first with a 
dummy interaction pulse with the same parameters as 
the real one, except for the value of the plateau current 
which should be slightly lower than the one of the real 
interaction pulses /^^^ such that no switching occurs in 
the interval Tp. Thus, the experimentalist can use this 
pulse to adjust the shape of the external flux pulse so 
that only switching events are recorded. This is helped as 
well by the fact that the changes Se{t) are second-order in 
the error representing the mismatch between the desired 
current pulse shape I{t) — —EjqcPq^ sin (/)(t) and the bias 
current signal that effectively reaches the sample, a direct 
consequence of the very definition of the optimal working 
point at which we are attempting to operate (first order 
changes in the qubit energy due to phase errors and fluc- 
tuations are zero). Also, the errors 5e{t) occur mostly 
during the relatively short raise and fall times (50 ns) 
and much less at the plateau - where the flux and the 
current being constant it is easier to ensure precisely the 
compensation. 

Finally, one can realize interaction-free experiments 
with even less control over the accumulated phase dif- 
ference A, using interaction pulses with slightly higher 
/(p) which tend to nullify predominantly the interfero- 
metric term. The price to pay for this is a reduction in 



the success rate. As an example, if I'^p^(I>o/Ejo = 96% 
and for a maximum error A = ±7r, the interference term 
becomes 7 times smaller than the interaction-free term 
(now also reduced to 0.27%). 

In conclusion, two fundamental physical processes, in- 
terferometry and tunneling, can be combined to demon- 
strate the equivalent of the non-classical interaction-free 
detection scheme for a superconducting quantum circuit. 
The crossover between standard interference effects and 
the interaction-free phenomenon is also discussed. 

This work was supported by the Academy of Finland 
(Acad. Res. Fellowship no. 00857 and projects no. 
7111994 and 7205476) and EU SQUBIT-2 (IST-1999- 
10673). The author is grateful to P. ZoUer for supporting 
a research visit at IQOQI Innsbruck. 



* Electronic address: |paraoanu@phys.jyTi.fi| 

[1] M. O. Scully and M. Suhail Zubairy, Quantum Optics 
(Cambridge University Press, Cambridge 1997). 

[2] Y. Makhlin, G. Schon, and A. Shnirman, Nature 386, 
305 (1999); J. E. Mooij et. al, Science 285, 1036 
(1999); Yu. Makhlin, G. Schon, and A. Shnirman, 
Rev. Mod. Phys. 73, 357 (2001). 

[3] Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Nature, 
398, 786 (1999); C. H. van der Wal et. al, Science 290, 
773 (2000); J. R. Friedman et. al, Nature 406, 43 (2000). 

[4] Yu. A. Paskin et. al, Nature 421, 823 (2003); T. Ya- 
mamoto et. al, Nature 425, 941 (2003); T. Duty, 
D. Gunnarsson, K. Bladh, and P. Delsing, Phys. Rev. 
B 69, 140503(R) (2004). 

[5] J. M. Martinis, S. Nam, J. Aumentado, and C. Urbina, 
Phys. Rev. Lett. 89, 117901 (2002). 

[6] I. Chiorescu, Y. Nakamura, C. J. P. M. Harmans, and 
J. E. Mooij, Science 299, 1869 (2003). 

[7] D. Vion at. al. Science 296, 886 (2002); A. Cottet et. al, 
Physica C 367, 197 (2002); D. Vion et. al Phys. Scr. T 
102, 162 (2002); A. Cottet, Ph. D. thesis, Univ. Paris 
VI (2002); E. Collin et. al, Phys. Rev. Lett. 93, 157005 
(2004); G. Ithier, et. al Phys. Rev. B 72, 134519 (2005). 

[8] A. Widom, J. Low Temp. Phys. 37, 449 (1979); 
A. J. Leggett and A. Garg, Phys. Rev. Lett. 54, 857 
(1985); A. J. Leggett, J. Phys.: Condens. Matter 14, 
R415 (2002). 

[9] A. Elitzur and L. Vaidman, Found. Phys. 23, 987 (1993); 
L. Vaidman, Found. Phys. 33, 491 (2003); L. Vaid- 
man, Quantum Opt. 6, 119 (1994); P. Kwiat et. al, 
Phys. Rev. Lett. 74, 4763 (1995); P. G. Kwiat et. al 
Phys. Rev. Lett. 83, 4725 (1999). 

[10] A. G. White, J. A. MitcheU, O. Nairz, and P. G. Kwiat, 
Phys. Rev. A 58, 605 (1998); J.-S. Jang, Phys. Rev. A 59 
2322 (1999); G. Mitchison and S. Massar, Phys. Rev. A 
63, 032105 (2001). 

[11] O. Hosten et. al Nature 439, 949 (2006). 

[12] J. Clarke et. al Science 239, 992 (1988); G. S. Paraoanu, 
Phys. Rev. B 72, 134528 (2005). 



